Whole-Genome Sequencing Analyses Reveal the Evolution Mechanisms of Typical Biological Features of Decapterus maruadsi.

Decapterus maruadsi is a typical representative of small pelagic fish characterized by fast growth rate, small body size, and high fecundity. It is a high-quality marine commercial fish with high nutritional value. However, the underlying genetics and genomics research focused on D. maruadsi is not comprehensive. Herein, a high-quality chromosome-level genome of a male D. maruadsi was assembled. The assembled genome length was 716.13 Mb with contig N50 of 19.70 Mb. Notably, we successfully anchored 95.73% contig sequences into 23 chromosomes with a total length of 685.54 Mb and a scaffold N50 of 30.77 Mb. A total of 22,716 protein-coding genes, 274.90 Mb repeat sequences, and 10,060 ncRNAs were predicted, among which 22,037 (97%) genes were successfully functionally annotated. The comparative genome analysis identified 459 unique, 73 expanded, and 52 contracted gene families. Moreover, 2804 genes were identified as candidates for positive selection, of which some that were related to the growth and development of bone, muscle, cardioid, and ovaries, such as some members of the TGF-ß superfamily, were likely involved in the evolution of typical biological features in D. maruadsi. The study provides an accurate and complete chromosome-level reference genome for further genetic conservation, genomic-assisted breeding, and adaptive evolution research for D. maruadsi.

Liver transcriptome analysis reveal the metabolic and apoptotic responses of Trachinotus ovatus under acute cold stress.

Trachinotus ovatus is an economically important fish and has been recommended as a high-quality aquaculture fish breed for the high-quality development of sea ranches in the South China Sea. However, T. ovatus shows intolerance to low temperature, greatly limiting the extension of farming scale, reducing production efficiency in winter, and increasing farming risks. In this study, liver transcriptome analysis was investigated in T. ovatus under acute low temperature conditions (20 and 15 °C) using RNA sequencing (RNA-Seq) technology. Inter-groups differential expression analysis and trend analysis screened 1219 DEGs and four significant profiles (profiles 0, 3, 4, and 7), respectively. GO enrichment analysis showed that these DEGs were mainly related to metabolic process and cell growth and death process. KEGG enrichment analysis found that DEGs were mainly associated with lipid metabolism, carbohydrate metabolism, and cell growth and death, such as gluconeogenesis, glycolysis, fatty acid oxidation, cholesterol biosynthesis, p53 signaling pathway, cell cycle arrest, and apoptotic cell death. Moreover, protein-protein interaction networks identified two hub genes (FOS and JUNB) and some important genes related to metabolic process and cell growth and death process, that corresponding to enrichment analysis. Overall, gluconeogenesis, lipid mobilization, and fatty acid oxidation in metabolic process and cell cycle arrest and apoptotic cell death in cell growth and death process were enhanced, while glycolysis, liver glycogen synthesis and cholesterol biosynthesis in metabolic process were inhibited. The enhancement or attenuatment of metabolic process and cell growth and death process is conducive to maintain energy balance, normal fluidity of cell membrane, normal physiological functions of liver cell, enhancing the tolerance of T. ovatus to cold stress. These results suggested that metabolic process and cell growth and death process play important roles in response to acute cold stress in the liver of T. ovatus. Gene expreesion level analysis showed that acute cold stress at 15 °C was identified as a critical temperature point for T. ovatus in term of cellular metabolism alteration and apoptosis inducement, and rewarming intervention should be timely implemented above 15 °C. Our study can provide theoretical support for breeding cold-tolerant cultivars of T. ovatus, which is contributed to high-quality productions fish production.

DataColor: unveiling biological data relationships through distinctive color mapping.

In the era of rapid advancements in high-throughput omics technologies, the visualization of diverse data types with varying orders of magnitude presents a pressing challenge. To bridge this gap, we introduce DataColor, an all-encompassing software solution meticulously crafted to address this challenge. Our aim is to empower users with the ability to handle a wide array of data types through an assortment of tools, while simultaneously streamlining parameter selection for rapid insights and detailed enhancements. DataColor stands as a robust toolkit, encompassing 23 distinct tools coupled with over 600 parameters. The defining characteristic of this toolkit is its adept utilization of the color spectrum, allowing for the representation of data spanning diverse types and magnitudes. Through the integration of advanced algorithms encompassing data clustering, normalization, squarified layouts, and customizable parameters, DataColor unveils an abundance of insights that lay hidden within the intricate relationships embedded in the data. Whether you find yourself navigating the analysis of expansive datasets or embarking on the quest to visualize intricate patterns, DataColor stands as the comprehensive and potent solution. We extend the availability of DataColor to all users at no cost, accessible through the following link: https://github.com/frankgenome/DataColor.

Whole-Genome Sequencing Analyses Reveal the Whip-like Tail Formation, Innate Immune Evolution, and DNA Repair Mechanisms of Eupleurogrammus muticus.

Smallhead hairtail (Eupleurogrammus muticus) is an important marine economic fish distributed along the northern Indian Ocean and the northwest Pacific coast; however, little is known about the mechanism of its genetic evolution. This study generated the first genome assembly of E. muticus at the chromosomal level using a combination of PacBio SMRT, Illumina Nova-Seq, and Hi-C technologies. The final assembled genome size was 709.27 Mb, with a contig N50 of 25.07 Mb, GC content of 40.81%, heterozygosity rate of 1.18%, and repetitive sequence rate of 35.43%. E. muticus genome contained 21,949 protein-coding genes (97.92% of the genes were functionally annotated) and 24 chromosomes. There were 143 expansion gene families, 708 contraction gene families, and 4888 positively selected genes in the genome. Based on the comparative genomic analyses, we screened several candidate genes and pathways related to whip-like tail formation, innate immunity, and DNA repair in E. muticus. These findings preliminarily reveal some molecular evolutionary mechanisms of E. muticus at the genomic level and provide important reference genomic data for the genetic studies of other trichiurids.

Integrated Transcriptomics and Metabolomics Reveal Changes in Cell Homeostasis and Energy Metabolism in Trachinotus ovatus in Response to Acute Hypoxic Stress.

Trachinotus ovatus is an economically important mariculture fish, and hypoxia has become a critical threat to this hypoxia-sensitive species. However, the molecular adaptation mechanism of T. ovatus liver to hypoxia remains unclear. In this study, we investigated the effects of acute hypoxic stress (1.5 ± 0.1 mg·L-1 for 6 h) and re-oxygenation (5.8 ± 0.3 mg·L-1 for 12 h) in T. ovatus liver at both the transcriptomic and metabolic levels to elucidate hypoxia adaptation mechanism. Integrated transcriptomics and metabolomics analyses identified 36 genes and seven metabolites as key molecules that were highly related to signal transduction, cell growth and death, carbohydrate metabolism, amino acid metabolism, and lipid metabolism, and all played key roles in hypoxia adaptation. Of these, the hub genes FOS and JUN were pivotal hypoxia adaptation biomarkers for regulating cell growth and death. During hypoxia, up-regulation of GADD45B and CDKN1A genes induced cell cycle arrest. Enhancing intrinsic and extrinsic pathways in combination with glutathione metabolism triggered apoptosis; meanwhile, anti-apoptosis mechanism was activated after hypoxia. Expression of genes related to glycolysis, gluconeogenesis, amino acid metabolism, fat mobilization, and fatty acid biosynthesis were up-regulated after acute hypoxic stress, promoting energy supply. After re-oxygenation for 12 h, continuous apoptosis favored cellular function and tissue repair. Shifting from anaerobic metabolism (glycolysis) during hypoxia to aerobic metabolism (fatty acid ß-oxidation and TCA cycle) after re-oxygenation was an important energy metabolism adaptation mechanism. Hypoxia 6 h was a critical period for metabolism alteration and cellular homeostasis, and re-oxygenation intervention should be implemented in a timely way. This study thoroughly examined the molecular response mechanism of T. ovatus under acute hypoxic stress, which contributes to the molecular breeding of hypoxia-tolerant cultivars.

OmicsSuite: a customized and pipelined suite for analysis and visualization of multi-omics big data.

With the advancements in high-throughput sequencing technologies such as Illumina, PacBio, and 10X Genomics platforms, and gas/liquid chromatography-mass spectrometry, large volumes of biological data in multiple formats can now be obtained through multi-omics analysis. Bioinformatics is constantly evolving and seeking breakthroughs to solve multi-omics problems; however, it is challenging for most experimental biologists to analyse data using command-line interfaces, coding, and scripting. Based on experience with multi-omics, we have developed OmicsSuite, a desktop suite that comprehensively integrates statistics and multi-omics analysis and visualization. The suite has 175 sub-applications in 12 categories, including Sequence, Statistics, Algorithm, Genomics, Transcriptomics, Enrichment, Proteomics, Metabolomics, Clinical, Microorganism, Single Cell, and Table Operation. We created the user interface with Sequence View, Table View, and intelligent components based on JavaFX and the popular Shiny framework. The multi-omics analysis functions were developed based on BioJava and 300+ packages provided by the R CRAN and Bioconductor communities, and it encompasses over 3000 adjustable parameter interfaces. OmicsSuite can directly read multi-omics raw data in FastA, FastQ, Mutation Annotation Format, mzML, Matrix, and HDF5 formats, and the programs emphasize data transfer directions and pipeline analysis functions. OmicsSuite can produce pre-publication images and tables, allowing users to focus on biological aspects. OmicsSuite offers multi-omics step-by-step workflows that can be easily applied to horticultural plant breeding and molecular mechanism studies in plants. It enables researchers to freely explore the molecular information contained in multi-omics big data (Source: https://github.com/OmicsSuite/, Website: https://omicssuite.github.io, v1.3.9).

Chromosome-Level Genome Assembly Provides Insights into the Evolution of the Special Morphology and Behaviour of Lepturacanthus savala.

Savalani hairtail Lepturacanthus savala is a widely distributed fish along the Indo-Western Pacific coast, and contributes substantially to trichiurid fishery resources worldwide. In this study, the first chromosome-level genome assembly of L. savala was obtained by PacBio SMRT-Seq, Illumina HiSeq, and Hi-C technologies. The final assembled L. savala genome was 790.02 Mb with contig N50 and scaffold N50 values of 19.01 Mb and 32.77 Mb, respectively. The assembled sequences were anchored to 24 chromosomes by using Hi-C data. Combined with RNA sequencing data, 23,625 protein-coding genes were predicted, of which 96.0% were successfully annotated. In total, 67 gene family expansions and 93 gene family contractions were detected in the L. savala genome. Additionally, 1825 positively selected genes were identified. Based on a comparative genomic analysis, we screened a number of candidate genes associated with the specific morphology, behaviour-related immune system, and DNA repair mechanisms in L. savala. Our results preliminarily revealed mechanisms underlying the special morphological and behavioural characteristics of L. savala from a genomic perspective. Furthermore, this study provides valuable reference data for subsequent molecular ecology studies of L. savala and whole-genome analyses of other trichiurid fishes.

Hiplot: a comprehensive and easy-to-use web service for boosting publication-ready biomedical data visualization.

Complex biomedical data generated during clinical, omics and mechanism-based experiments have increasingly been exploited through cloud- and visualization-based data mining techniques. However, the scientific community still lacks an easy-to-use web service for the comprehensive visualization of biomedical data, particularly high-quality and publication-ready graphics that allow easy scaling and updatability according to user demands. Therefore, we propose a community-driven modern web service, Hiplot (https://hiplot.org), with concise and top-quality data visualization applications for the life sciences and biomedical fields. This web service permits users to conveniently and interactively complete a few specialized visualization tasks that previously could only be conducted by senior bioinformatics or biostatistics researchers. It covers most of the daily demands of biomedical researchers with its equipped 240+ biomedical data visualization functions, involving basic statistics, multi-omics, regression, clustering, dimensional reduction, meta-analysis, survival analysis, risk modelling, etc. Moreover, to improve the efficiency in use and development of plugins, we introduced some core advantages on the client-/server-side of the website, such as spreadsheet-based data importing, cross-platform command-line controller (Hctl), multi-user plumber workers, JavaScript Object Notation-based plugin system, easy data/parameters, results and errors reproduction and real-time updates mode. Meanwhile, using demo/real data sets and benchmark tests, we explored statistical parameters, cancer genomic landscapes, disease risk factors and the performance of website based on selected native plugins. The statistics of visits and user numbers could further reflect the potential impact of this web service on relevant fields. Thus, researchers devoted to life and data sciences would benefit from this emerging and free web service.

Gene Expression Profile and Co-Expression Network of Pearl Gentian Grouper under Cold Stress by Integrating Illumina and PacBio Sequences.

Pearl gentian grouper (Epinephelus fuscoguttatus ♀ × Epinephelus lanceolatus ♂) is a fish of high commercial value in the aquaculture industry in Asia. However, this hybrid fish is not cold-tolerant, and its molecular regulation mechanism underlying cold stress remains largely elusive. This study thus investigated the liver transcriptomic responses of pearl gentian grouper by comparing the gene expression of cold stress groups (20, 15, 12, and 12 °C for 6 h) with that of control group (25 °C) using PacBio SMRT-Seq and Illumina RNA-Seq technologies. In SMRT-Seq analysis, a total of 11,033 full-length transcripts were generated and used as reference sequences for further RNA-Seq analysis. In RNA-Seq analysis, 3271 differentially expressed genes (DEGs), two low-temperature specific modules (tan and blue modules), and two significantly expressed gene sets (profiles 0 and 19) were screened by differential expression analysis, weighted gene co-expression networks analysis (WGCNA), and short time-series expression miner (STEM), respectively. The intersection of the above analyses further revealed some key genes, such as PCK, ALDOB, FBP, G6pC, CPT1A, PPARα, SOCS3, PPP1CC, CYP2J, HMGCR, CDKN1B, and GADD45Bc. These genes were significantly enriched in carbohydrate metabolism, lipid metabolism, signal transduction, and endocrine system pathways. All these pathways were linked to biological functions relevant to cold adaptation, such as energy metabolism, stress-induced cell membrane changes, and transduction of stress signals. Taken together, our study explores an overall and complex regulation network of the functional genes in the liver of pearl gentian grouper, which could benefit the species in preventing damage caused by cold stress.

Complete mitochondrial genome of yellowback seabream, Dentex hypselosomus and phylogenetic analysis of the family Sparidae.

The complete mitogenome sequence of Dentex hypselosomus was amplified by designing 15 primer pairs. The circle genome was 16,618 bp in size and the overall base composition was 27.04% of A, 26.38% of T, 17.08% of G, and 29.50 of C, with significant anti-G bias. The complete mitogenome of D. hypselosomus encodes 37 canonical mitochondrial genes, two non-coding regions, an L-strand replication region (OL), and a control region (D-loop). The D-loop contained termination sequence domain (TAS), central conserved domains (CSB-F, CSB-E, CSB-D, CSB-C, and CSB-A), and conserved sequence blocks (CSB-1, CSB-2, and CSB-3). Phylogenetic analysis of nine sparid species well supported the phylogenetic position of D. hypselosomus and revealed the phylogenetic relationship of the family Sparidae at the level of mitochondrial genomes.

Complete mitochondrial genome and phylogenetic relationship of ornate threadfin bream, Nemipterus hexodon (Perciformes, Nemipteridae).

The complete mitochondrial genome of the ornate threadfin bream, Nemipterus hexodon, was first determined by the pairs-walking sequencing in this study. The circular mtDNA molecule was 17,115 bp in size and the overall nucleotide composition of H-stand was A (29.55%), T (27.36%), G (16.08%), and C (27.01%), with a slight bias towards A + T. The complete mitogenome encoded 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 2 non-coding regions (an origin of L-strand replication and a control region). The Bayesian tree supported the phylogenetic position of N. hexodon, which provided useful information for phylogenetic relationship in genus Nemipterus.

The second complete mitogenome of Nemipterus virgatus to dissect control region structure and phylogenetic problem of the superfamily Sparoidea (Teleostei, Perciformes).

In this study, we determined the complete mitogenome of Nemipterus virgatus of which the length was 17,073 bp, including 37 canonical mitochondrial genes and 2 non-coding regions. The control region contained termination associated sequence domain (TAS), central conserved domain (CSB-F, CSB-E, CSB-D, CSB-C, and CSB-A), conserved sequence block domain (CSB-1, CSB-2, and CSB-3), and tandem repeat sequence domain (TTD). Nine single nucleotide polymorphisms and three insertion of tandem repeat sequence (each length in 28 bp) were detected between two N. virgatus mitogenomes. The phylogenetic analysis showed that the families Nemipteridae, Sparidae, Centracanthidae, and Lethrinidae did not gather into a monophyly of superfamily Sparoidea in the neighbor-joining tree.

The second complete mitogenome of Scolopsis ciliata (Perciformes, Nemipteridae) to analyze control region structure and intraspecific variation.

In this study, we reported the complete mitogenome of Scolopsis ciliata sampled from the Beibu Gulf, South China Sea. It was 16,734 bp in length and contained 37 canonical mitochondrial genes and 2 non-coding regions. Fifty-one single nucleotide polymorphisms and five insertion and deletions were detected between two S. ciliata mitogenomes, suggesting a high level of intraspecific genetic variation in the species. The control region contained termination associated sequence domain (TAS), central conserved domain (CSB-F, CSB-E, CSB-D, CSB-C, and CSB-A), and conserved sequence block domain (CSB-1, CSB-2, and CSB-3). The phylogenetic analysis of 21 complete mitogenome sequences well supports the phylogenetic position of S. ciliata and reveals the phylogenetic relationship of the genus Scolopsis.